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The decays B — > Kvv arise from flavor-changing neu- 
tral currents (FCNC), which are forbidden at tree level 
in the Standard Model (SM). The lowest-order SM pro- 
cesses contributing to these decays are the W box and 
the Z penguin diagrams shown in Fig. [TJ New physics 
contributions may enter at the same order as the SM. 
These contributions, some of which could increase the 
branching fraction by up to ten times relative to the SM, 
include: unparticle models Minimal Supersymmet- 
ric extension of the Standard Model at large tan ft [|J , 
models with a single universal extra dimension [3| , scalar 
Weakly Interacting Massive Particle (WIMP) dark mat- 
ter [l] and WIMP-less dark matter @. A recent SM 
prediction (ABSW model 0) for the total B -> Kvv 
branching fraction is (4.5 ± 0.7) x 10~ 6 , while an ear- 
lier prediction (BHI model [7j), based on a different form 
factor model, is (3.81^) x 10" 6 . The BHI model was 
used by previous analyses [8j, |9[ and provides a baseline 
for comparison between results. The current experimen- 
tal upper limit (UL) on the total branching fraction for 
B + —} K + vv (charge conjugation is implied throughout) 
is 1.4 x 10~ 5 at the 90% confidence level (CL) from the 
Belle Collaboration [8j, while an earlier BABAR analysis 
set an UL of 5.2 x 10" 5 (90% CL) @. The only existing 
UL on the total branching fraction for B° — » K°vv is 
1.6 x 10~ 4 (90% CL) from Belle 0. 

We report results of a search for B + — > K + vv and 
B° — > K°vv, with branching fractions for both decays 
as well as for the combination B — >■ Kvv. We also re- 
port on partial branching fractions for B + — > K + vv in 
two regions of di- neutrino invariant mass squared (q 2 )- 
The low-g 2 region (q 2 < QAm 2 B ) is selected by requiring 
p* K+ > 1.5GeV/c and the high-q 2 region (q 2 > 0Am 2 B ) 
by P*k+ < 1.5GeV/c in the T(AS) center-of-mass system 
(CMS) 10], where tub is the mass of the B meson and 
p* K+ is the magnitude of the CMS 3-momentum of the 
signal K + candidate. The high-g 2 region is of theoreti- 
cal interest because the partial branching fraction in this 
region could be enhanced under some new physics mod- 
els 0. 

This analysis is based on a data sample of (459.0 ± 
5.1) x 10 6 BB pairs, corresponding to an integrated lu- 
minosity of ~ 418 fb -1 of e + e~ colliding-beam data and 
recorded at the T(AS) resonance with the BABAR detec- 
tor [HI at the SLAC PEP-II asymmetric-energy B Fac- 
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FIG. 1: Lowest-order Feynman diagrams for B — ¥ Kvv, with 
the W box on the left and Z penguin on the right. 



tory. Charged particle tracking is provided by a five-layer 
silicon vertex tracker and a 40-layer drift chamber in a 1.5 
T magnetic field. A CsI(Tl) electromagnetic calorimeter 
(EMC) is used to measure photon energies and directions 
and to identify electrons. All quantities in this paper 
which are measured by the EMC are required to exceed a 
minimum 20 MeV cluster energy, unless a higher thresh- 
old is explicitly noted. The magnetic flux return from 
the solenoid, instrumented with resistive plate chambers 
and limited streamer tubes (IFR), provides muon identi- 
fication. We identify K + candidates by using a detector 
of internally reflected Cherenkov light (DIRC) as well as 
ionization energy loss information from the tracking sys- 
tem. 

Due to the presence of two neutrinos in the B — > Kvv 
final state, it is not possible to exploit the kinematic 
constraints on the B mass and energy which are typi- 
cally used to distinguish signal and background events 
in B meson decays at the T(AS). Instead, before look- 
ing for the signal decay, we first reconstruct a B decay 
(B Tec ) in one of several exclusive D^*Hv semileptonic final 
states. We then search for the signal B — > Kvv among 
the remaining charged and neutral particles in the de- 
tector that are not part of the B rcc . We collectively re- 
fer to these remaining particles as B Toe for "rest of the 
event." This strategy is common to several BABAR anal- 
yses [lH EH and has the advantage of higher efficiency 
compared with reconstruction of the B rec in hadronic de- 
cay modes 0- 

We reconstruct the D candidates in the following de- 
cay modes: K~ir + , K~tt + it + , K~7r + ir + ir~ , K~ir + 7T°, 
K®ir + , and K^tt + tt^ . The K® candidates, reconstructed 
in the K® — > tt + tt~ mode, are required to have a ir + n~ 
invariant mass within 25 MeV/c 2 of the nominal K® mass. 
D candidates are similarly required to have a recon- 
structed invariant mass within 60 MeV/c 2 of the nominal 
value [l4j, except for the K ~n + ir° mode where the range 
is 100MeV/c 2 . We form D*° D°tt°, D*+ D°tt+, 
and D* + — > D + ir° candidates with a required mass dif- 
ference (m(D*) - m{D)) in the range 130-170 MeV/c 2 . 
In addition, we combine D and 7 candidates to form D* 
candidates with a required mass difference in the range 
120-170 MeV/c 2 . A DM candidate is combined with 
an identified electron or muon with momentum above 
0.8 GeV/c in the CMS to form a B lec candidate. In 
events with multiple reconstructed B rQC candidates, we 
select the candidate with the highest probability that the 
daughter tracks originate from a common vertex. After a 
B rcc candidate has been identified, the remaining charged 
and neutral decay products are used to classify the -B r oe 
as either a background event or a possible signal candi- 
date. 

As a first step in refining the selection of B me can- 
didates, we veto K candidates which, when combined 
with a remaining charged or neutral pion candidate, 
have a Kw invariant mass within 75 MeV/ c 2 of the nom- 
inal if* (892) mass. We also veto events where a re- 
maining charged track can be combined with a ir° can- 
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didate to yield a p + candidate, with a mass window 
0.45 < m(p + ) < 1.10GeV/c 2 . Similarly vetoed are 
events where three remaining charged tracks can be com- 
bined to yield an candidate, with a mass window 
0.6 < m(a~i) < 2.0GeV/c 2 . These vetoes eliminate, with 
little loss of signal efficiency, sizable backgrounds that 
consist mostly of random track combinations. After the 
vetoes, B + (B°) signal candidate events are required to 
possess K + (K® —> tt + it~) candidates, accompanied by 
at most two (one) additional charged tracks, which are 
assumed to have been incorrectly left out of the B rcc . For 
the K + final state, the B lec lepton daughter and the K + 
are also required to be oppositely charged. For the Kg 
final state, signal candidates are required to have a tt + tt~ 
invariant mass within 25 MeV/c 2 of the nominal K® mass. 

At this stage of the selection, each event has a S rec 
candidate representing a B meson reconstructed in a 
semileptonic decay and a B roc candidate formed from 
the rest of the event, with the latter representing the sig- 
nal decay. In simulated K + (Kg) signal events that have 
passed this selection, 99% (92%) of events have a cor- 
rectly identified signal K + (Kg). However, a large back- 
ground still remains. Further background suppression is 
achieved using a multivariate event selection algorithm, 
a bagged decision tree (BDT) [TH, [l6[ , that can leverage 
many weak discriminating variables to achieve high back- 
ground rejection. Such an algorithm needs to be trained 
with simulated signal and background events, henceforth 
referred to as Monte Carlo (MC) events. We use a 
GEANT4 [13 detector simulation to obtain large samples 
of simulated signal events generated with a pure phase- 
space model (which are later re-scaled to the BHI signal 
model), as well as samples of non-resonant e + e~ — > qq 
(q = u, d, s, c), BB, and t + t~ background events, whose 
sizes are one (uds), two (cc), three (BB), and one (t + t~) 
times luminosity. These background events are aug- 
mented with a separate sample, with a size 13 times lu- 
minosity, of simulated BB doubly semileptonic events, 
the largest source of background. 

We construct two ensembles of BDTs, one for the K + 
signal mode and one for the K®. To create an ensem- 
ble, we repeatedly divide the total signal and background 
datasets in half randomly, creating 20 distinct BDT 
training and validation datasets, where each dataset has 
a 50% correlation with any other because approximately 
50% of the events are shared. This procedure makes op- 
timal use of the limited statistics of MC events that pass 
the initial event selection and results in a more statisti- 
cally precise unbiased estimate of background contribu- 
tions. Use of the ensemble of 20 BDTs created for each 
final state also averages out the variations in BDT re- 
sponse compared to a single BDT trained and validated 
with a single division of the simulated signal and back- 
ground datasets [HI, [l9|. The choice of 20 divisions, in- 
stead of a lower or higher number, represents a balance 
between minimizing the variation versus minimizing the 
overhead of multiple BDTs. 

Each BDT of the K+ (K°) ensemble uses 26 (38) dis- 



criminating variables, described in the Appendix. These 
variables fall into four general categories: quantities re- 
lated to the missing energy in the event, to the overall 
event properties, to the signal kinematics, and to the 
overall reconstruction quality of the B Tec . Some quan- 
tities are given in two different frames and thus allow 
the BDTs to extract from them additional discriminat- 
ing power. Several additional variables were initially con- 
sidered but were pruned during the BDT optimization 
process, as they were found to add little additional sen- 
sitivity. 

"Missing Energy" quantities relate to the fact that sig- 
nal events are expected to possess significant missing en- 
ergy and momentum because the signal decay includes 
two neutrinos. In contrast, the dominant background 
events usually acquire missing energy and momentum 
as a result of particles passing outside of the detector 
fiducial acceptance, with the result that distributions of 
quantities related to missing energy differ between signal 
and background. 

After the B ICC and K or K® signal candidate have 
been identified, signal events are expected to have lit- 
tle or no additional activity in the detector, other than a 
few low-energy clusters in the calorimeter resulting from 
hadronic shower remnants, beam backgrounds, or simi- 
lar sources. In contrast, background events arising from 
higher-multiplicity B decays typically possess additional 
charged or neutral particles within the detector. Vari- 
ables which characterize this additional detector activ- 
ity can provide discriminating power between signal and 
background, and are indicated by the term "extra" in the 
following. 

The strongest discriminant for both K + and K a en- 
sembles is _E extra , the sum of all detector activity not 
explicitly associated with either the B TCC or K signal can- 
didate, followed by p* K+ for the K + ensemble and by the 
lab energy of the signal K® for the K° ensemble. The 
reconstructed mass of the D from the J5 rec is the third 
ranking variable for both channels. 

Figure shows signal, background, and data distribu- 
tions from the validation set of K + and K° BDT output 
for a BDT randomly selected from the 20 BDTs in the 
ensemble. The other 19 BDTs are similar to that shown. 

We choose as the target signal efficiency the one that 
maximizes expected signal significance averaged over the 
20 BDTs, under the assumption of a branching fraction 
of 3.8 x 10~ 6 . This signal significance is s/Vs + b, where 
s is the number of signal events, and b is the number 
of background events. Optimization using a figure of 
merit based upon signal efficiency and independent of 
assumed branching fraction yields similar results. For 
each BDT, a BDT output value that yields the target 
signal efficiency is calculated. For example, the BDT 
output cuts for the BDTs shown in Figure [5] are 0.976 
for the K+ BDT and 0.955 for the K° BDT. The mean 
background for target signal efficiency is obtained by av- 
eraging the individual background estimates from each of 
the 20 BDTs. Thus, we treat each ensemble of 20 BDTs 
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FIG. 2: (a) K + and (b) K° BDT output for data (diamonds), background MC (solid), and signal MC (dotted) events. For 
each plot, the scale for the data and background events is on the left axis, and the scale for the signal events is on the right 
axis. The distribution of signal MC events is normalized to unit area. 



as a set of correlated estimators for the numbers of signal 
and background events in a signal region defined by the 
target signal efficiency. 

The low-q 2 (high-q 2 ) measurement uses the K + en- 
semble but only includes events with p* K+ > 1.5GeV/c 
{p*k+ < 1.5GeV/c), which means that only those events 
are used to calculate the signal efficiency and the back- 
ground prediction. The low-g 2 measurement has the 
same BDT output cuts and background prediction as the 
primary K + measurement, with only the signal efficiency 
changed by the restriction on p* K+ . On the other hand, 
the high-g 2 measurement has its own set of BDT output 
cuts based upon its own optimized signal efficiency, along 
with its own background prediction. 

The total optimized signal efficiency for the K + (K°) 
mode is 0.16% (0.06%), while the efficiency for the K + 
low-g 2 (high-g 2 ) region is 0.24% (0.28%). The uncer- 
tainty in the signal efficiency is discussed below. Figure 
[3] shows the BDT selection efficiency versus p* K for the 
K + ,K°, and high-g 2 measurements, where the BDT se- 
lection efficiency considers only the effect of the BDT 
output cut. 

To measure the branching fractions, we use the value 
obtained from simulated events of the predicted back- 
ground in the signal region, the number of observed data 
events, and the signal efficiency, as shown by the fol- 
lowing equation: B = (N b s — Nbk g )/e.NB, where B is 
the branching fraction, N a i, s is the number of observed 
data events, N^kg is the number of predicted background 
events, e is the total signal efficiency, and Nb is the num- 
ber of B mesons, either charged or neutral [2fj| . that 
are relevant to the branching fraction. We account for 
the 50% correlation between each of the datasets when 
computing the statistical uncertainty of the estimated 



background contribution by using a standard method for 
combining correlated uncertainties [l9j . 

Data control samples are used to ensure that both 
signal-like and background-like events in actual data are 
classified similarly to simulated events. The vetoed af 
events offer a high-statistics control sample which can 
be used to compare the K + and K° BDT output dis- 
tributions for background events in both simulated and 
actual data. We find good agreement between data and 
MC events in the BDT output distribution for both final 
states, with only a (+5± 2)% data-MC discrepancy. For 
the K + mode we make a +5% adjustment to the expected 
number of background events, based upon a weighting 
technique that corrects data-MC discrepancy in the side- 
band K + BDT output next to the signal region, and we 
assign the full adjustment as a systematic uncertainty. 
Likewise, for the high-g 2 K + measurement, we make a 
+25% correction to the expected number of background 
events and assign the full correction as a systematic un- 
certainty. In the K® final state, we find a (+10 ± 3)% 
data-MC discrepancy in the sideband BDT output next 
to the signal region, and we make a +10% correction and 
assign the full correction as a systematic uncertainty. 

To validate our signal efficiency estimates and assess 
their systematic uncertainties, we use high-purity sam- 
ples of B+ -> K+J/iP{-> £+£-) decays (where i+£~ = 
e + e~, fi + [i~). The two leptons from the J/ip are dis- 
carded in order to model the unseen neutrinos of the 
signal decay, and then the events are subjected to the 
same selection requirements as other signal candidates. 
Classifying J/ipK data and MC events, we find only 
a (-10 ± 10)% data-MC discrepancy in the BDT out- 
put distribution. Although we do not make any cor- 
rection, we assign a 10% systematic uncertainty to the 
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FIG. 3: BDT selection efficiency in the signal region versus p* K for (a) K + , (b) Jfg , and (c) high-g 2 K + simulated signal events, 
considering only the effect of the BDT output cut. 



TABLE I: Systematic uncertainties 



Category 



Uncertainty 



Signal efficiency 
K + background prediction 
High-g 2 K + background prediction 
Ks background prediction 



14% 
5% 
25% 
10% 



TABLE II: Total signal efficiencies and MC expectations of 
the number of data events. The uncertainties shown are sys- 
tematic for iVsgni, with statistical negligible, and statistical 
followed by systematic for iVbkgd- 



Mode 




e (in %) 


iVsgnl 


-/Vbkgd 


K+ 




0.16 


2.9 ±0.4 


17.6 ±2.6 ±0.9 


K° s 




0.06 


0.5 ±0.1 


3.9 ± 1.3 ±0.4 


low-g 2 


K+ 


0.24 


2.9 ±0.4 


17.6 ± 2.6 ± 0.9 


high-g 


2 K+ 


0.28 


2.1 ±0.3 


187 ± 10 ± 46 



estimated signal efficiency for all four measurements 
(K + ,K®, low-q 2 K + , high-q 2 K + ) based on these results. 
We also assign a signal efficiency systematic uncertainty 
of 10% to account for the theoretical uncertainties of the 
signal models. Adding these in quadrature, we assign a 
total uncertainty of 14% in the estimation of signal effi- 
ciency for both final states. Table U summarizes all of the 
systematic uncertainties. 

Table |TT] shows the total signal efficiencies and the ex- 
pected number of signal and background events in the 
data. We performed a blind analysis where data events 
with BDT outputs above the optimized values were not 
counted or plotted until the analysis methodology and 
sources of systematic uncertainty were fixed as described 
above. 

Table Hill shows our results. The non-integer number of 
observed events results from averaging the integer yields 
from the 20 BDTs of each type. We calculate two-sided 
68% confidence intervals for the number of excess events 
based on the statistical and systematic uncertainties in 
the background estimates, as well as the statistical un- 
certainty on the number of events observed in the data. 



TABLE III: Observed and excess data events, with statistical 
uncertainties [2lJ shown for iVobs and combined statistical and 
systematic uncertainties shown for N exceaa . The last column 
shows the probability that excess events could be due solely 
to a background fluctuation. 



Mode 




N obs 


N 

1 * excess 


Prob. 


K+ 
K° 

\ow-q 2 
high-g 


K+ 
2 K + 


19.4+lt 
6.1+*'° 

mill 


9+4.1 
Z - Z -2.S 

1 s+6.2 

90+49 
zo -48 


38% 
23% 
38% 
33% 



TABLE IV: Branching fraction (BF) central values and upper 
limits. The low- and high-g 2 values are partial BFs, while the 
rest are total BFs. 



Mode 




BF 


90% CL 95% CL 






xlO -5 


xlO -5 


xlO -5 


K+ 
K° 




2 +u - 8 
u.z_ 7 

1 7 +3.i 


1.3 
5.6 


1.6 
6.7 


Comb. 


K+,K° 


n 5+ - 7 

u.o_ 7 


1.4 


1.7 


Low-g 2 


K+ 


2 +0 - 6 

U ' Z -0.5 


0.9 


1.1 


High-g 2 


K+ 


1 C+3-8 
_i - 8 -3.8 


3.1 


4.6 



Figure [4] shows the averaged BDT outputs in the signal 
region for K + , K°, and high-g 2 K + data overlaid with 
the background and signal contributions, while Figure [5] 
shows similar plots for the p* K distribution in the signal 
region. Figure [B] shows the integrated numbers of events 
(observed, predicted background, and excess over back- 
ground) in the signal region for K + , K°, and high-g 2 
K+ data for each of the 20 BDTs of each type. Table [TV] 
gives the branching fraction central values, along with 
corresponding 90% and 95% CL upper limits, assuming 
the BHI signal model (the ABSW model gives similar re- 
sults) . The upper limits are calculated using a frequentist 



method [21]. The quoted uncertainties include all statis- 
tical and systematic uncertainties. Our results constrain 
the B — > Kvv branching fraction at the 90% CL to a few 
times the SM expectation, with limits of 1.3 x 10~ 5 for 
B+ -> K+vV and 5.6 x 10~ 5 for B° -> K°vV. 

We are grateful for the extraordinary contributions of 
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FIG. 4: Averaged BDT signal-region output for (a) K + , (b) K%, and (c) high-q 2 K + data, with expected signal and background 
contributions. The signal estimate assumes a branching fraction of 3.8 x 10 -6 . 
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FIG. 5: Averaged p* K signal-region output for (a) K + , (b) K® , and (c) high-q 2 K + data, with expected signal and background 
contributions. The signal estimate assumes a branching fraction of 3.8 x 10 -6 . 
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FIG. 6: Integrated numbers of observed (red triangles), expected background (black circles), and excess events (blue squares) 
for data for each BDT: (a) K + ', (b) Kg, and (c) high-q 2 K + . The individual uncertainties are purely statistical and assume no 
correlation between data sets. The horizontal dashed lines show the sum of the statistical and systematic uncertainties on the 
mean number of excess events. 
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Appendix 

Definitions of BDT Variables 

In the following the notation [K + ] or [K°] indicates that 
a variable is used only by that ensemble; otherwise it is 
used by both BDT ensembles. 

BDT input variables related to missing 4-momcntum 
The event missing 4-momentum is computed from the 
difference between the 4-momentum of the combined 
e + e~ beams and the 4-momenta of all charged and neu- 
tral particles reconstructed in the detector. 

• Energy component of missing momentum 4-vector 

• Energy component of missing momentum 4-vector 
(CMS) 

• Magnitude of the missing momentum 3- vector 

• Magnitude of the missing momentum 3-vector 
(CMS) 

• Cosine of the angle with respect to the beam axis 
of the missing momentum 3-vector 

• Cosine of the angle with respect to the beam axis 
of the 3-momcntum vector representing the differ- 
ence between the initial event momentum and the 
summed momenta of the B rcc and B loc candidates 
[K°] 

BDT input variables related to overall event properties 

• -E"extra = ZJiEi, where Ei is the energy of an isolated 
EMC cluster or a charged track and the sum is over 
all tracks or clusters which are not part of the B lec 
or the B Ioe 

• Total energy of all reconstructed charged and neu- 
tral particles in the event 

• Minimum invariant mass obtained from the combi- 
nation of any three charged tracks in the event 

• Total charge of all tracks in the event [K°] 



• Total charge of all tracks matched to EMC energy 
deposits [K°] 

• Number of extra EMC clusters 

• Number of Kl candidates in the EMC 

• Number of IFR K L candidates [K+] 

• Number of extra reconstructed tracks 

• Magnitude of the 3-momentum of a candi- 
date T(AS) computed from the B rcc and B ioc 
4-momenta [K°] 

• Angle with respect to the beam axis of a candidate 
T(45) 3-momentum vector computed from the B rec 
and B roe 4-momenta [K°] 

• Normalized second Fox- Wolfram moment of the 
overall event 

BDT input variables related to signal kinematics 

• Cosine of the angle between the signal K and the 
event thrust axis 

• Cosine of the angle between the signal K and the 
Dl thrust axis 

• Energy of the signal kaon [K°] 

• Reconstructed invariant mass of the signal K® [K°] 

• Magnitude of the 3-momentum of the signal kaon 

• Magnitude of the CMS 3-momentum of the signal 
kaon 

• Cosine of the angle with respect to the beam axis 
of the 3-momentum vector of the signal kaon 

• Uncertainty in the x-component of the signal K 
point of closest approach to the e + e~ interaction 
point, as determined from a three dimensional fit, 
with the x-axis defined perpendicular to the beam 
axis in the horizontal plane of the detector [K°] 

• Uncertainty in the cc-component of the signal K 
point of closest approach to the e + e~ interaction 
point, as determined by a fit in the xy-plane, 
with the x-axis defined perpendicular to the beam 
axis (z) in the horizontal plane of the detector [K°] 

BDT input variables related to B mc reconstruction 

• x 2 P er degree of freedom of the vertex fit of the 
tracks making up the B rcc 

• co sg by = (2E* eam ■ E* m - m|„ - m 2 m ) / (2 P * Dl ■ 

V E tea m ~ m Bn)' where E beam * One half the to- 
tal CMS energy, niBn is the nominal B meson 
mass [HI and E* Dl , itidi and p* Dl are the CMS en- 
ergy, invariant mass and 3-momentum magnitude 
of the D - lepton combination used in the recon- 
struction of the -B rec 
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• cos by re-calculated with the addition of a photon 
to the Dlv candidate such that 100 < (m(D ,^) — 
m{D )) < 150MeV/c 2 

• Reconstructed decay mode of the D from the B rcc 

• Uncertainty in the x-component of the point of clos- 
est approach to the e + e~ interaction point of the 
leading pion daughter from the D meson, with the 
x-axis defined perpendicular to the beam axis in 
the horizontal plane of the detector 

• Number of daughters possessed by the recon- 
structed D from the B Tec 

• Number of extra 7r° candidates satisfying 0.115 < 
m(77) < 0.150 GeV/c 2 and E 1 > 30 MeV 

• Reconstructed invariant mass of the B lcc 



• Reconstructed invariant mass of the D from the 

-^rcc 

• Magnitude of the CMS 3-momentum of the B rec 

• Magnitude of the CMS 3-momentum of the D can- 
didate from the B rec 

• Magnitude of the 3-momentum of the lepton from 
the B vcc [K°] 

• Magnitude of the CMS 3-momentum of the lepton 
from the B TCC 
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